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Abstract. We report on the theoretical investigation of photonic crystal cavities 
etched on a suspended membrane for the generation of polarization entangled photon 
pairs using the biexciton cascade in a single quantum dot. The implementation of 
spontaneous emission enhancement effect increases the entanglement visibility, while 
the concomitant preferential funncling of the emission in the cavity mode increases the 
collection of both entangled photons. We demonstrate and quantify that standard 
cavity designs present a polarization dependent emission diagram, detrimental to 
entanglement. The optimization of HI cavities allows to obtain both high collection 
efficiencies and polarization independent emission, while keeping high Purcell factors 
necessary for high quality entangled photon sources. 
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Entangled photon sources play an important role in quantum communication 
networks or quantum information processing |ll[2]- For the former, they are a necessary 
resource for the realization of quantum repeaters [3] based on quantum teleportation or 
quantum entanglement swapping. In the first demonstrations of such relays, parametric 
down conversion sources have been used for the generation of entangled photon pairs 
[H El [6l [7]. Such non-linear sources of entanglement can combine narrow spectral 
bandwidths with a maximal generation rate [3 [H [9]. Although these sources may 
be very useful and easy to implement, they always suffer from the Poissonian statistics 
of the emitted photons pairs leading to multipair emission and thus decreasing the 
fidelity of entanglement [lOj. Being able to produce polarization entangled photon pairs 
on demand would be an important step towards robust quantum relays. Such sources 
can be obtained from the biexciton-exciton cascade emission of a single quantum dot 
[TT] . and first experimental demonstrations have been reported [T2l [T3] . Obtaining 
entangled photons pairs, however, from such quantum dots sources with both high 
fidelity and high collection efficiency remains a problem. Implementing Cavity Quantum 
ElectroDynamics effects by embedding a single quantum dot in a microcavity could not 
only improve the fidelity of the emitted pair |,14J by taking advantage of the Purcell 
effect, but also by enhance the collection efficiency [16] . One promising microcavity 
for such purpose is the single defect hole cavity in a triangular lattice of holes (HI) 
etched on a suspended membrane, due to its small mode volume and its polarization 
degeneracy. However, in a standard HI cavity, the radiation pattern of the two 
fundamental degenerate modes do not overlap, leading to photon distinguishability and 
thus destroying entanglement. Theoretical calculations demonstrate that this radiation 
pattern can be strongly modified by changing, for instance, the position of the holes 
surrounding the defect [17]. This modification of the design is necessary to avoid 
distinct emission patterns. This emission pattern distinguishability is related to the 
mode overlap. 

In this paper, we report on the theoretical investigation of HI photonic crystal 
cavities etched on slab membrane, in order to obtain both high collection efficiencies 
for both photons and a high overlap between the two fundamental energy-degenerate 
modes. The dependency of the Bell inequalities as a function of the mode overlap is 
derived. We also investigate the impact of the position of the quantum dot inside the 
cavity on entanglement visibility and collection efficiency. 

1. Entangled state density matrix for non overlapping modes 

Polarization entangled photon pairs can only be obtained if and only if, even in principle, 
the polarization of the photon can not be determined by measuring another degree of 
freedom as for example the photon's energy. In the same way, if the emission mode 
of one of the photons of the pair does not perfectly match the emission mode of the 
other photon, the non-maximal overlap between the two emission modes will reduce the 
fidelity of entanglement. Our analytical derivation of this non- maximal mode overlap 
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effect is based on the density matrix of tlie plioton pair emitted by tlie cascade emission 
from tlie biexcitonic level of a single quantum dot. 



Energy 




Figure 1. Schematic description of the two-photon cascade in a typical quantum 
dot four-level system with an energy splitting 2hduj of the relay level, yielding two 
collinearly polarized photons (either H or V). 



The eigenbasis of the dot involves four levels : |2) (biexcitonic level), \1h) and 
(two excitonic levels with opposite angular momenta) and |0) (fundamental level). In 
this eigenbasis, the emitted photons are linearly polarised along the horizontal (H) or 
vertical (V) directions. The density matrix of the photon pair in this particular basis 
B = [H1H2, H1V2, V1H2, V1V2] where the subscript z = 1, 2 is related to the photon 
emitted by the biexcitonic level and excitonic level respectively, can be written in the 
form |I14|: 
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with 7i the exciton decay rate. F = Yu -|- Fy is the cross-dephasing rate between the 
two excitonic states. F j^jp ± 5F describe phenomenologically relaxation mechanisms 
between the two excitation states, leading to incoherent population transfers between 
these two states {8T takes into account the possible rate asymmetry in this process). 
buj is the energy splitting of the excitonic levels. Fi is the pure dephasing rate induced by 
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dephasing processes that occur simultaneously and attach the same information to the 
phase and energy of these two excitonic levels. and Ty are polarization-dependent 
pure dephasing rates induced by dephasing processes that do not identically affect the 
two relay levels and whose impact depends on the polarization of the excitonic states. 

We consider that the exciton and biexciton photons are emitted in the same 
cavity mode, ie the cavity mode is resonant with both transitions. The cavity mode 
is doubly-degenerate in polarisation, due to the Cg symmetry of the HI cavity. We can 
consequently describe the photons polarization in the {H, V) basis defined previously, 
independently of the orientation of the dot with regards to the orientation of the photonic 
crystal. We define $_ff(r) (resp. $y(r)) the complex spatial far field distribution of the 
horizontal [H) (resp. vertical {V)) polarization modes. Propagation occurs along the 
orthogonal direction to the photonic crystal membrane and f denotes the radial vector 
perpendicular to the propagation axis. The density matrix can be rewritten under the 
form: 

pin, r~'2)xy,uv = "^xinY * ^yir~'2y * Pxy,uv * ^ui^) * $^(r"^2) (3) 

with xy and uv E B and 'Pxy,uv being the density matrix element on line xy and column 
uv. p{fi,f2)xy,uv is the density matrix element on line xy and column uv of the new 
density matrix ^(ri,?^). Let t(r) be the function describing the detectors' active areas 
which are placed along the propagation axis. There are in fact two distinctive detectors 
(one for each photon of the pair [T3]) but we suppose that they have the same sensitive 
area for the sake of simplicity. The density matrix can be reduced for the detected 
photon pairs to: 

^ _ /(iVi (iVa t{rl)t{r~*2)p{rl,f*2) 
Tr(/ (i^n d'^r2 t(rl)t(r^)p(rl, r^)) 
We consider that the cavity is positioned at the focal point of a microscope objective, 
which transforms the emitted far field into the complex transverse shape of a propagative 
beam. In the first order approximation, and $y are real and positive, corresponding 
to the case where the transverse phase is constant in the propagative modes (plane wave 
approximation). Let k and e be: 



k= I d^r t(r)^$j^(f)$v(r) (5) 
'd\ t{r)^jj{r) = [ d^r t{r}^l{r} (6) 



The overlap factor K can be expressed as K = k'^/e'^ and the final expression of the 
detected photon density matrix in the case of non-maximal overlap of the two emission 
modes is: 

fa {d-ici)K\ 

i - a C2K 

C2K I - a 

V {d + ici)K a J 

Note that only the coherence terms are modified, and are multiplied by the overlap 
factor K. When both modes do not overlap {K = 0), the mutual coherence is 
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(7) 



Optimizing HI cavities for the generation of entangled photon pairs 



5 



erased and entanglement vanishes. On the contrary, maximally entangled states can 
only be obtained for K = 1. Following [H], Bell inequalities can be rewritten as 
S = 2^2 (a + K*{d- cs)) > 2. Even in the case of a single dot emitting maximally 
entangled photons, a minimum overlap of K > 2/\/2 — 1 = 41% is required in order to 
violate Bell inequalities. 

2. HI cavity for rricLximally entangled photons 

We consider the HI cavity as potential candidate for the generation of entangled photon 
pairs since it sustains two energy degenerate dipole modes with a field maximum in the 
center of the cavity. This cavity offers both a low mode volume and theoretically high 
quality factors by fine tuning the inner holes |T7j . 




(a) (b) 

Figure 2. HI photonic crystal cavity. The inner holes are shifted by the quantity d 
from (standard HI design) to 0.18a. The blue line is the plane where the field is 
registered at the end of the simulation for radiative pattern calculations. 

Simulations of HI cavities were performed with the 3D finite-difference time-domain 
(FDTD) method, using a freely available software package with subpixel smoothing 
for increased accuracy[TS]. The simulated structure is depicted in figure [21 The HI 
photonic crystal (figure [2](a)) has a lattice constant equal to a =270 nm and the holes 
have a radius of r^ =80 nm. The refractive index of the GaAs membrane is equal 
to n =3.46. Above and below the membrane, a free space volume is inserted, with 
a thickness of 3a. The simulation volume is finally surrounded by Split Field Perfect 
Matched Layers (PML). A temporally short Gaussian dipole pulse (with a width of 10 
optical oscillations) is launched in the center of the cavity (figure [2]|^b)) and used as a 
white light source. After extinction of the source, the electromagnetic field evolves freely 
over a time corresponding to approximately 300 optical cycles of the source, after which 
all low quality factor modes have radiated, thus leaving only the desired cavity mode 
in the simulation volume. In such conditions, the decay of the field amplitude at some 
fixed non-nodal point inside the cavity follows a simple exponential function of rate Fg. 
The emission wavelength is determined by measuring the optical oscillation frequency. 

The collection efficiency is defined as the ratio between the incident power within a 
given emission cone normal to the membrane, over the total emitted power. At any given 
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time, the total emitted power is given by P^e/ = T' * W where W denotes the energy 
inside the cavity, and T' is the intensity decay rate {V = 2rc). Let U be the total energy 
in the simulated volume at time t. The energy outside the cavity (Uout) corresponds to 
the energy emitted by the cavity which has not yet reached the edge of the simulation 
volume, Uout = V *W * D/c with D being the radius of the simulated volume and c 
the speed of light. Thus U = W * {1 + D/c). Since D is only a few micro-meters wide, 
D/c « 1, and the total emitted power can be written as Pref = SFcf/. 

The emission mode of the cavity is estimated following reference [19] (mainly Eq. 
23). This method relies on the complex value of the electromagnetic field on a plane 
(P) just above the membrane (see figure [2](b)) at some time (in our case the end of the 
simulation). The real part of the electromagnetic field is directly measured on (P), and 
the imaginary part is deduced from measurements of the real part of the electromagnetic 
field a quarter oscillation later (at the cavity's resonant frequency), taking into account 
the losses induced during this quarter of cycle. This allows us to extract the far field 
emission mode from the light-cone of the spatial Fourier transform of the field with a 
unique simulation run in real values, thus saving valuable calculation time. Emission 
patterns for a membrane thickness oik = 0.26/xm as a function of the hole displacement 
d are depicted in figure [3] and correspond qualitatively with the emission patterns 
calculated by Roemer et al [T7j using a 3D Finite Element Maxwell Solver. For small 
hole displacements {d < 0.10), the emission diagram is almost spherical, whereas by 
increasing the hole displacement {d > 0.14), a pronounced, directional Gaussian-like, 
central peak appears. 



d=0 d=0.09a d=0.12a 




Figure 3. Emission patterns for various hole displacement d for a membrane thickness 
of h =0.26/im. Each pattern is normalized to its maximum value. The distance 
to the center is sin{9) where 6 is the normal angle to the membrane. The gray 
(resp. black) circle represents an objective of numerical aperture NA=0.5 (resp. 
NA=sm(7r/4) = 0.7) 

A systematic analysis of the maximal Purcell factor {F^°'^), collection efficiency 
(rj), and mode overlap {K) has been performed by varying two parameters: the hole 
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Maximal Purcell factor Collection efficiency Mode overlap 
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Figure 4. Maximal Purcell factor, collection efRciency, and mode overlap as a function 
of the hole displacement (d, in cristal units a) and membrane thickness {h, in /im). 



displacement d and the membrane thickness h. Figure H] summarizes the results where 
the mode overlap and collection efficiency are estimated for an objective with a numerical 
aperture of 0.5. By maximal Purcell factor, we mean the Purcell factor obtained for a 
resonant punctual monochromatic source placed at the maximum of the electrical field 
intensity. A more detailed view of the variation of these parameters as a function of 
the hole displacement is given in figure [5l in which the membrane thickness is fixed to 
h = 0.26/im and three numerical apertures have been taken into account (0.2, 0.5 and 
0.7). The wavelength Ac of the cavity depends strongly on both parameters d and h. A 
good linear approximation is Xc = d x 0.28 fim + h x 0.69 + 0.82/im with 4 nm of maximal 
deviation. A global homothetic transformation of the design, including the membrane 
thickness, should be latter applied to tune the cavity to the desired wavelength but is 
not taken into account here. 

As we demonstrated earlier [T3] a Purcell factor of 10 should be sufficient to restore 
entanglement in the emission cascade from single quantum dot. From this point of 
view alone, the whole domain of variation of the two parameters studied here satisfies 
this condition. However, spontaneous emission enhancement of the exciton rate is not 
sufficient to realize a deterministic efficient entangled photon source: efficient coupling 
of both photons to the cavity mode, high emission mode overlap and high collection 
efficiency are also needed. Let us ffist consider the coupling of both photons to the 
cavity mode. Exciton and biexciton lines are usually separated by about 2nm. If 
the exciton is resonant with the cavity in order to obtain the highest spontaneous 
emission acceleration, the biexciton should preferably be also quasi-resonant with the 
cavity mode, in order to beneffi from the redirection of the emission and to increase the 
collection efficiency of the biexcitonic photon. If we set a minimum Purcell factor for the 
quasi- resonant biexcitonic line to a value of about 5, this in turn limits the Purcell factor 
which can be reached for the excitonic line to a maximum value of about 110 for a cavity 
of modal volume 0.7(Ac/^)'^. This limitation excludes a large domain of parameters 
value around the hole shift of d = 0.13. Let us consider now the problem of emission 
mode overlap and collection efficiency. Since the highest collection efficiencies with a 
Gaussian-like emission pattern are reached for a membrane thickness of h = 0.26/im, 
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Hole displacement (a) Hole displacement (a) Hole displacement (a) 

Figure 5. Maximal Purcell factor, collection efficiency, and mode overlap as a function 
of the hole displacement (d, in cristal units a) for a membrane thickness h = 0.26^to. 
The collection efficiency and the mode overlap are calculated for numerical apertures 
NA=0.2 (red curves), NA=0.5 (green curves) and NA=0.7 (blue curves). 

At low hole displacement [d < 0.09a), high mode overlap can be obtained (in 
excess of 95%) but at the expense of a low collection efficiency (below 10% for standard 
microscope objectives (NA< 0.5) ) and does not exceed 22% for high numerical aperture 
objectives (NA=0.7) with a theoretical maximum of 50% since light is emitted upwards 
and downwards with the same intensity. The emission diagram (figure [3]) for a hole 
displacement of d=0.1a clearly explains the situation. The mode is almost uniform in 
every direction giving rise to a high mode overlap, and the collection efficiency scales 
as the objective's numerical aperture. In the range 0.11a < d < 0.14a the mode 
overlap presents a distinct dip down to 57% for an objective with NA=0.7, and the 
collection efficiency drops to the same extent. In this region, the Purcell effect reaches 
its maximum (d=0.145a). The drop in the collection efficiency is correlated to the 
increase of the quality factor, corresponding to a better confinement of the light inside 
the photonic crystal slab and a reduction of vertical losses at the F point. For larger 
hole displacements {d > 0.15a) the collection efficiency increases sharply reaching 22% 
for a NA=0.5 (and d = 0.16a), corresponding to a 4-fold increase compared to a hole 
displacement of = 0.10a. At the same time, the mode overlap increases up to 96% 
reaching almost the values obtained at low values of d. The mode profile (fig [3]) is 
almost TEqo in the propagation direction perpendicular to the membrane. A numerical 
aperture of NA=0.2 increases the overlap up to almost 100% but at the expense of a 
low collection efficiency. On the other hand, a numerical aperture NA=0.7 increases 
the collection efficiency by a factor of 1.5 compared to a numerical aperture NA=0.5, 
but the mode overlap does not exceed 83% indicating that almost half of the energy is 
astigmatic. As a conclusion, a hole displacement of d = 0.16a with an objective with 
numerical aperture of NA=0.5 and membrane thickness h = 0.82Ac/n seems to be the 
optimum in terms of collection efficiency (of the order of 22% with a maximal calculated 
value of 50%) and mode overlap (of the order of 96%). 
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3. Impact of the position of the dot in the cavity 

Until now the quantum dot has been considered to be perfectly placed in the center 
of the cavity, implying that both polarizations undergo the same Purcell effect and 
that the cavity mode is equally fed for both polarizations. Deterministically aligning 
a photonic crystal around a single quantum dot so that the dot is positioned in the 
center of the cavity is technologically extremely challenging but mandatory. Several 
techniques are being developed [Sni EI] although due to experimental uncertainties 
which are essentially due to the electronic beam lithography process, the mismatch 
of the quantum dot position with respect to the center of the HI cavity can be up to 10 
nm. 

The position mismatch implies a breaking of the Cq symmetry. The position of 
the dot will be identified by a direction X, as shown on figure Et^a), the Y direction 
being orthogonal to the X direction. The two polarization modes of the cavity remain 
unchanged. Therefore the sustained modes of the cavity will be described in the basis 
{X,Y) and no more in the {H,V) basis. Due to the mismatch with the cavity modes, 
the dipole will preferentially excite one of the modes {X or Y polarized) leading to an 
unbalance of the fraction j3i [i G [X, F]) of spontaneous emission in the cavity mode. 
Inevitably, this will in turn impact the entanglement visibility. 




(a) (b) (c) 

Figure 6. (a): Normalized /? factor as a function of the position mismatch of the 
dot. The red (blue) curve corresponds to the X (Y) polarization, (b): Figure of merit 
(r = T{"'"^5a;/(?iF™"^)) as a function of the position mismatch for various values of 
the Bell test: S=2, S=2.2, S==2.4, S=2.6 and S=2.8. 

Figure [6](b) depicts the normalized jSi factor as a function of the position mismatch 
along the X direction for both polarizations. The fraction of spontaneous emission for 
both polarizations {f3x and /3y) at null mismatch; this permits us to normalize the /3j 
to 1 when the dot is centered, and use in our simulations the extrapolated amplitude of 
the corresponding mode divided by the amplitude of the whole field at the maximum of 
the excitation, in order to deduce the evolution of its 

The effect of an asymmetric Purcell factor on the Bell's inequality is modeled in 
the annexe of this paper. If we want to reach some value 5" of the Bell test whereas 
the dot is misplaced, this puts a maximal limit to a figure of merit defined as the 
adimensional ratio r = Tl'^^^dui / [hF^'^^) where T^"^^^ is the bulk lifetime of the dot, 
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2hSLV the excitonic energy splitting and F^°-^ the maximal Purcell effect of the cavity 
(at zero mismatch). Entanglement visibility increases when r tends to zero. This figure 
of merit only depends on the shortest lifetime that can be obtained in the dot-cavity 
system and on the quantum beat period between the two exciton states. Figure [6] (c) 
depicts the evolution of this figure of merit for different S values, calculated using the 
model introduced in the Annexe. For example, in the case of a centered dot with an 
excitonic bulk lifetime of 1 ns, an excitonic energy splitting of about 2fieV, submitted 
to a maximal Purcell factor of 10 and not subjected to incoherent processes, the figure 
of merit r = 0.3 allows S to reach a value above 2.6. Conversely, the Bell's inequality 
is hardly violated (5* = 2) if the same dot is about 70 nm away from the center of the 
cavity. For a more usual value of the excitonic splitting (bfieV), the maximal mismatch 
enabling for Bell's inequality violation drops to 10 nm. 

4. Conclusion 

In this paper we derived the Bell inequalities for a quantum dot in a photonic crystal 
cavity with non-overlapping polarization modes. By analyzing the emission pattern of 
modified HI cavities, we demonstrate that it is possible to obtain both high collection 
efficiencies (of the order of 22%) and maximally overlapping modes while keeping high 
Purcell factors. Finally we estimate that the position of the quantum dot with respect to 
the cavity center has to be more accurate than 50 nm, in order to implement an efficient 
quantum dot source of polarization entangled photons from the biexciton cascade. 
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Annexe: effect of an assymetry in the Purcell effect 

The asymmetric branching ratio induced by a polarisation dependant Purcell factor can 
be modelled as follows. Let the state of the system (dot and optical fields) be 

|^(t)) =P2(t)|2;0;0)+ ^ / c/cu2P„(o;2, t)|l„; m^, cua; 0) 

u=H,V 

+ J2 duj2dujipuu{(^i,(^2,t)\(];lt,u;2;lt,uji) (8) 

u=H,v' 

where the first of the three entries within the ket refers to the quantum dot's level, the 
two other entries refer to the first and second emitted photons of polarisation Tt and 
pulsation uji (i=l, 2 respectively) (see figure [H We distinguish here the emission rates 71 
and 72 with respect to both polarizations H and V. We assume that incoherent processes 
are negligible, so that crossed terms combining horizontal and vertical orientations 
disappear.) 
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The expressions of the p2, Pu and Puu coefficients are determined using the Wigner- 
Weisskopf approximation. Considering the system at long times (t >> I/71 and 72) the 
terms P2(t) and pu{uj2,t) tend to zero, which gives a state that can be factorized into 
a radiative part I'^r) and the fundamental source state |0). Considering no spectral 
filtering, the density matrix of the photon pair in the polarization basis is: 

duJiduJ2{uJ2,UJi\\'$R){'$R\\uJ2,UJi) (9) 



E 

u,v=Hy 



1 



itlt) {if if 



2(1 + 2(5F2) 



/ {SF + lf 





duJiduJiPuui^i, ^2, oo)p™(t^i, ^2, 00)* 



l+ig 










l-ig 




i6F-lf ) 



(10) 



(11) 



where we defined the relative difference of Purcell factors 5F = {Fh — Fy)/{FH + Fy) 
and the normalized splitting g = 25uj / {'-^^^^{Fh + Fy)). SF is approximated from the 
ration between the modal coupling factors as: {Ph — Pv)/{.Ph + Pv)- In the same way 
as in [H], we deduce the expression S of the Bell test, from which we deduce data 
presented on figure [6](c). 
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